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The memory kernel of the self-part of the intermediate scattering function is studied. We found 
that the short-time behavior of the memory kernel for our MD-simulated glass- forming Ni2oZr8o- 
system shows a deviation from the polynomial form predicted by MCT. By using the gaussian 
approximation we give a more detailed description of the short-time dynamics of the system. 
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I. INTRODUCTION 

At present the phenomena behind the liquid-glass tran- 
sition and the nature of the glassy state are not fully un- 
derstood, despite the progresses of recent years. In con- 
trast to usual phase transformations the glass transition 
seems to be primarily dynamic in origin and therefore 
new theoretical approaches have to be developed for its 
description. Several theoretical models have been pro- 
posed to explain the transition and the corresponding 
experimental data. The latter concern of both the tem- 
perature dependence of particular properties, such as the 
shear viscosity and the structural relaxation time, and 
the time dependence repsonse of the spectra as visible in 
the dielectric susceptibility, inelastic neutron scattering, 
and light scattering spectra investigations. The spectral 
measurements have been extended to cover the large fre- 
quency range from below the primary a-relaxation peak 
up to the high-frequency region of microscopic dynamics 
dominated by vibrational modes [l[ . 

One of the promising approaches in this field is mode 
coupling theory (MCT). MCT originally was developed 
to model critical phenomena 0, Q . The non-classical be- 
havior of the transport properties near the critical point 
was thought to be caused by nonlinear couplings between 
slow (hydrodynamics and order parameter) modes of the 
system. In later years, the MCT was found to be applica- 
ble more generally, to describe nonlinear effects in dense 
liquids 0] and nonhydrodynamic effects in the case of the 
glass transition. 

The idealized schematic model of the MCT for the 
liquid-glass transition, firstly proposed by Bengtzelius 
et.al. [5| and independently by Leutheusser [Qj, implies 
that the liquid (" ergodic" ) state is characterized through 
a decay of structural correlations while below a critical 
temperature Tc a solidlike ("nonergodic") state exists 
with structural arrest where fluctuations are effectively 
frozen. In this state correlations decay on a longer time 
scale by thermally activated diffusion processes 0, [a, Q . 

According to the MCT the time dependence of the cor- 
relation functions for structural fluctuations follows from 
a damped harmonic oscillator equation with a retarted 



damping term. All the details concerning the time evo- 
lution of the structural correlations are hidden in the 
memory kernel that describes the retarted damping. 

In a previous paper [l3|, hereafter refered to as I, one 
of us has evaluated this memory kernel from the time de- 
pendence of the correlation functions in the molecular dy- 
namics(MD) simulations of the glass-forming NisgZrso- 
system. Here two results from I shall be recalled: first, it 
found that the short-time behavior of the memory kernel 
of the self-part of the intermediate scattering function 
<&((/, t) shows a deviation from the polynomial form com- 
monly assumed in the idealized (or extended) schematic 
model of the MCT. This deviation, in particular, explains 
the absence of the inverse powerlaw decay of the corre- 
lation predicted by the MCT. The deviations originate 
from the vibrational motions of the atoms in local poten- 
tial minima and are restricted to the short times, below 
some 0.3 ps. This, on the other hand, explains that the 
MCT predictions for the late /3-regime and for the a-peak 
are in a good agreement with the experimental results 
[H [13, H B lH and MD-simulation ^^. Second, 
there is a characteristic quantity (?,„, as introduced in 
[18J . that can be used as a relative measure to describe 
how close a liquid structure has approached the arrested 
nonergodic state. If gm > 1 for a given temperature T, 
then the system is able to undergo the transition into 
an idealized nonergodic state with a structural arrest; if 
gm < 1, the system remains in the (undercooled) liquid 
state (T > Tc), and if gm ~ 1, the system stands in a 
critical situation. 

The behavior of the memory kernel F{t) and the pa- 
rameter gm are formally related through the character- 
istic function g(<I>) = F{t){l/^{q,t) - 1) introduced in 
[18| . gm nieans the maximum of ^('I') for $ S [0, 1] and 
the short time behavior of F{t) is, in particular, visible in 
the behavior of g{^) for values <&((/, t) ^ 1. Therefore we 
here are interested in a detailed description of (?($) for 
short times, respectively $((?,t) — > 1. Since, according 
to I, <&((/, i) for short times is determined by the vibra- 
tion spectrum of the amorphous structure, the memory 
kernel F(t) for short times and the corresponding g also 
depend on the vibration spectrum. The latter is visi- 
ble in the velocity autocorrelation function (VACF) and 



thus we in the following look for an interrelationship be- 
tween the VACF and g. Moreover, we here consider the 
Nii_a;Zr2; system at a; = 0.8, while in I the concentration 
a; = 0.5 was investigated. Thus we ask to which extend 
the results for x — 0.8 confirm the results for x = 0.5 
from I. 

Our paper is organized as follows: In Section |TT1 we 
present the model and give some details of the com- 
putations. Section |TTT] gives a brief discussion of some 
aspects of the mode coupling theory as used here and 
it describes the gaussian approximation to relate $(5,^) 
and the VACF at the short time. Results of our MD- 
simulations and their analysis are then presented and 
discussed in Section [IVl 



II. SIMULATIONS 

As in paper I, the present simulations are carried out 
as state-of-the-art isothermal- isobaric ( A^, T, p) calcula- 
tions. The Newtonian equations of A^ = 648 atoms (130 
Ni and 518 Zr) are numerically integrated by a fifth or- 
der predictor-corrector algorithm with time step Ai = 
2.5xl0~^^s in a cubic volume with periodic boundary 
conditions and variable box length L. With regard to 
the electron theoretical description of the interatomic 
potentials in transition metal alloys by Hausleitner and 
Hafner [1^ we model the interatomic couplings as in [23| 
by a volume dependent electron-gas term E^oiiV) and 
pair potentials 4'{r) adapted to the equilibrium distance, 
depth width, and zero of the Hausleitner-Hafner poten- 
tials [131 for Ni2oZr8o [21|. For this model simulations 
were started through heating a starting configuration up 
to 2000 K, yielding a homogeneous liquid state. The 
system then is cooled continuously to various annealing 
temperatures with cooling rate —dT = 1.5x10^^ K/s. Af- 
terwards the obtained configurations of various annealing 
temperatures (here 1500-800 K) are relaxed by modelling 
an additional isothermal annealing that depend on the 
object temperature. Finally the time evolution of these 
relaxed configurations is modelled and analyzed. Full de- 
tails of the simulations are given in [21| . 



III. THEORY 

A. Mode Coupling Theory 

In the MCT equation of motion for the correlation 
function of structural fluctuations with wave vector q can 
be expressed as a damped harmonic oscillator equation 
(see, e.g, [Hil) 



$(g,t)/r!2 + <i>(g,t) 



dTF{q,t-T)^q,T) = Q (1) 



with initial conditions $(9,0) — 1 and $(g,0) — 0. fio is 
a microscopic (phonon) frequency and F{t) is the mem- 
ory kernel. The schematic modell assumes that the time 



dependence of F{q,t) can be expressed by that of ^q{t) 
and the idealized version of MCT, neglecting atomic dif- 
fusion, relies on the assumption 



F{q,t) = Fo{t) -.^ h{t) + fmq,t)) 



(2) 



where /($) means the polynomial in $ and h{t) a short 
time viscous damping which conveniently is approxi- 
mated by an instantaneous term r]QS(t). The asymptotic 
behavior of ^{q,t) is determined by /($). 

As mentioned above, there is for a given temperature 
a characteristic function g{^) whose maximum value can 
be used to determine whether the system is able to un- 
dergo a structural arrest ( gm > 1) or whether it re- 
mains in the (undercooled) liquid state ( gm < 1 )• This 
characteristic function can be derived as follows: Given 
lim $ w $°° > 0. Let us assume that there is a critical 

t — >CXD 

time tr so that for t > tr 



m = , 



$(i) = 0. 



(3) 



(4) 



Substitute the equations © and (g]) in MCT-Eq. H]), 
then the MCT-Eq. becomes 



$°° + /($°°)($°°_l) = 0, 
and after a little algebra we find 

.9($) = /(<D)(l/<D-l) 



(5) 



(6) 



with $ G [0,1]. This function is related in simple way 
to AF($) = /($) - <I)/(1 - $) frequently used in the 
schematic MCT [2^]. 

According to the extended version of MCT which sim- 
ulates atomic diffusion by taking into account the cou- 
pling to transverse currents, the memory kernel has the 
following form 



F{t)=C-'{l/[D{z) + 



^Fo{t)), 



-]} 



(7) 



Here C means the Laplace transform, C^^ its inverse. 
D{z) models the coupling to the transverse currents. 
F{t) from Eq.([7]) leads to the final decay of structural 
fluctuations also below Tc- This sets one of the basic 
problems in classification of the solutions of Eq. ([T]) as it 
may not be obvious for a given solution whether it be- 
longs to the regime above or below Tc if atomic diffusion 
is included. 

For our analysis we introduce the Laplace Transform 
for the memory kernel F{t) and the correlation function 
^(t), by taking z = e — iuj with e — > 0, as follows 



lim£{F(t)}, 

e— >0 

lim/:{$(t)k 



Fciiu) + iFsiiu) (8) 

$c(w) + z$,(w) (9) 



where $c('^) a-nd $s('^) are the cosines-part and sinus- 
part of the Fourier transformed $(t). It should be noted 



that according to the fluctuation-dissipation theorem the 
loss part of the dynamic susceptibiUty spectrum x"(w) is 
related to the dynamics structure factor $c(w) via 



X"(tj)=w$c(w) 



(10) 



Substitute equations ^ and Q in MCT-Eq. H]), then 
we obtain the equation for the spectral distribution 
uFc{uj) of the memory kernel 

coF^H = ^^ , . (11) 

[1-C^<1>,(C^)]V[C^$C(C^)]' 

Taking the inverse- Fourier transformed Fc{lu), we have 
the time dependent memory kernel 



Fit) 



du)Fc{Ljo) cos{uot) 



(12) 



B. Gaussian Approximation 

According to the theoretical models studied in the the- 
ory of hquids [2^, [2a, [23] we want to illustrate the correla- 
tion functions in the short-time regime for temperatures 
near the critical temperature. One of the models based 
on a cummulant expansion [28| express the correlation 
function in an isotropic system 



$((?,i) - exp {-q^piit) + q^p2{t) -...] 



where 



pi(i) = -{J|x,(i + to)-x,(to) 
n{t)f[R2{t)-l] , 

3 {J|x,;(i + to)-X,(io)| 



P2{t) = l[pi{t)f[R2it) 



R2{t) 



Xi(t +to) -xj(to)r 



(13) 

(14) 
(15) 

(16) 



The assumption that $(g, t ) ca n be expressed only by the 
leading term of expansion [2^, [23| 



$((7, t) = exp { -^ {\v{t + to)-r(to)|' 



(17) 



is called the Gaussian approximation. Also in this ap- 
proximation the correlation function is related to one- 
sixth the mean square displacement (MSD) of atoms dur- 
ing time t. In section IV. we will show our calculating 
results of Eq. (jl7p and prove that this approximation is 
correct in the short-time regime. 

To show, by using this Gaussian approximation, 
whether our system in the short-time regime is a purely 
harmonic system or not, we need a relation between the 
mean squared displacement and the spectral distribution 
Z{ijj) of the velocity autocorrelation function Zit) and a 



assumption that in purely harmonic system Z{lj) would 
be proportional to the density of state (DOS) Zdosii-^)- 
One find that relation 



|r(i + to)-r(to)r) -6w^ / dio 



(1 — C0s{L0t)) 



UJ^ 



where 



(18) 



'{uj) = / dtcos{ujt)Z{t) 
Jo 

{vi{t + to) • Wi(io)) 



dtcos{ujt)- 



Wj(io)r 



2 1 /l ^+ m2\ f^B 



M 



, (19) 



(20) 



Actually Ngai et.al.(see e.g. 29]) have developed a 
model, which is called " Coupling- Model" , in order to 
explain the dynamics of glass-forming liquids. Here we 
briefly review that model. The model made use of the 
assumption that vibration (phonon) and relaxation (dif- 
fusion) of the molecules contribute independently to the 
density-density correlation function; thus, $(g, t) is equal 
to the product ^phiq,t) x ^reiaxiq,t) [30] • In this model 
the correlation function of the relaxation ^reiax (q-, t) 



c{q,t) 



eM-t/to{q,T)) 
exp(~(t/T(g,r))i-") 



for t <tc 

for t > tn 



(21) 

where respectively tc is a temperature-independent 
crossover time separating two time regimes in which the 
relaxation dynamics differ, ig the Debye-relaxation time, 
T the Kohlrausch-relaxation time, and 1 — n — (3 the 
Kohlrausch stretching exponent 

The phonon contribution to $(g,t) is determined by 
the harmonic DOS of the phonon modes, Z^^™{lo), ac- 
cording to the formula (see e.g. [30, [sll), in a classical 
limit, 



$pft(g,i)=exp(-^w(i,T)), 



(22) 



where 



u{t,T) 



2kBT 
M 



'^-(i^^^Zr"(-). (23) 



u{t, T) decreases with time, leveling off to constant value 
u{T) = limu(i, r). Correspondingly, $(g,i) has the 

t — *oo 

value exp(— (7^u(r)/2), which is the well known Lamb- 
Mossbauer factor. For a single Einstein oscilator of fre- 
quency uje the harmonic DOS of the phonon modes is 
given by 



Ztr{u)^5{u:-u,) 
and, for a Debye model of a solid, 



yharm 
^doa 



iu) 



( Sut'^ /ujf) for U! < LlJ]J 

for uj < ujd 



(24) 



(25) 



By comparing this expressions of the couphng model 
with those of the theory of Uquids and making a assump- 
tion that Z^°J'"((jj) and Zdos{^) are the same quantities 
in the short-time regime we find that ^ph{q,t) fEq. iP^ 
with Eq. ((23l) ) of the couphng model is the same as $(g, t) 
(Eq.dni) with Eq.llHl)) of the theory of hquids (with the 
Gaussian approximation). 

Considering that the model of the theory of liquid has 
a simple form we use for our further analysis that model. 
We believe that by making use of the model it is enough 
to prove whether in the short-time regime our system is 
a harmonic system or not. 

Following Cummins et.al.[l| the correlation func- 
tions $(t) and the corresponding susceptibility spectrum 
x"(w) can be considered to consist of three regimes: 1) 
at high frequencies (short times) a nearly temperature- 
independent microscopic structure; 2) at low frequen- 
cies (long times) a strongly temperature-dependent a- 
relaxation peak associated with structural relaxation and 
diffusion process; 3) (intermediate frequencies) between 
a) and b) a minimum in x"(^) {ot "plateau") in $(t) 
whose amplitude and position is also temperature de- 
pendent. 

From above definition of regions we write now in gen- 
eral form the full susceptibility spectrum 

x"{T,uj) = x'st{T,oj) + x'!t{T,uj) + XitiT,^) (26) 

After our intensive study about the short-time regime for 
the object temperatures near the critical temperature it 
is a little difficult to separate purely the short time regime 
from the intermediate regime. Also it is important to 
keep in mind that our definition for the short-time part in 
Ea. (p6| include a small part of the /3- relaxation process. 
Considering that the spectral distribution of the veloc- 
ity autocorrelation function Z(t) for each temperature 
near T^ ( from [19] Tc = 1050 ± 25K ) shows too in 
our system a nearly temperature-independent form (see 
Fig. [T]), we can scale now the short-time part of the sus- 
ceptibility spectrum, so that 






(27) 



where x'^j is the master curve of the susceptibility spec- 
trum that is now a temperature-independent quantity 
and is defined 



1 ^ T 



(28) 



where N is the number of the object temperature and Tq 
is a arbitrary normalization of temperature. 
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FIG. 1: The spectral distribution of the velocity autocorrela- 
tion function for different temperature near T^. 



The bracket means averages over the atom i and the ini- 
tial configurations tg. 

In Figs. [2] and [3] we show the results evaluated from our 
MD data (with symbols) for wave vector q = ^(n,0, 0) 
with n = 8 and 6 which correspond respectively to j^l = 
19.7 nin~^ and \q\ = 14.7 nm~^. Both figures present the 
average over the star of six q vectors which are equivalent 
on assuming equivalence of the simulation cube axis and 
their inverses. 



A. MCT-Analysis 

The following Fig. [5] and [3] present the self-part of the 
intermediate scattering function for different object tenr- 
peratures according to MCT (with line) and our MD- 
simulation (with symbols). The discrepancy between 
both results are considerable as results from the bad fit- 
ting of the Kohlrausch-law to our MD-data and the inac- 
curracy by calculation (see [10|, [iSl, |2]| for details of the 
algorithm used for calculation). 
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IV. RESULTS AND DISCUSSION 

The central object of our analysis is the self-part of the 
intermediate scattering function: 



$(q,i) = (exp{iq- [x,(t + to) - Xj(io)]}) 



(29) 



FIG. 2: The self-part of the intermediate scattering function 
for wave vector qe, from MD-data (with symbol) and from 
MCT-Equations (with line). 

From both figures we found that the self-part of the in- 
termediate scattering function show for gg and gg at aver- 
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FIG. 3: The same as Fig. [2] but for qe. 



age temperatures the structural relaxations that happen 
in three succesive steps. According to predictions of the 
idealized schematic model of MCT the correlator should 
decay in three succesive steps [20] at this temperatures. 
The first is a fast initial decay on the time scale of vi- 
brations of atoms {t < 1 ps). This step is characterized 
by MCT only global. The second is the /3-relaxation 
regime (typically in the range 1 ps < i <1 ns). In the 
early /^-relaxation regime the correlator should decrease 
according to $(i) = fc+A/f^ and in the late /3- relaxation 
regime, which appear only in melting, according to von 
Schweidler-law fc — Bt^ . Between them the wide plateau 
appear near the critictical temperature Tc. In the case 
of idealized glass the von Schweidler-law does not ap- 
pear more, but <l>(i) = fc + A/f^ flows into the kon- 
stant asymptotic value fc- In melting the a- relaxation 
appear as the last step that results from the decay af- 
ter the von Schweidler-law and could be described by the 
Kohlrausch-law ^{t) = Aoexp(— (i/ro)^) whose the re- 
laxation time To near the glass transition shifts dramatic 
to the longer time scale. 

The decay for the early /3-regime, which according 
to MCT follows the inverse power-law ^ ^ fc + A/t"-, 
could not be observed in our data. The reason for it 
is that in our system $m(r) (the maximum of g(<&)) 
and $o(r) (the threshold-value of G($), see Eq.(12Sl)) are 
nearly parallel, which cause that the power-law of the 
early /3-regime is dressed up. On the one side, for this 
regime the so-called "Boson-peak" is discussed by sev- 
eral authors (see e.g.fs^, [23, [SJ), on the other side one 
called this peak as the /3-maximun [331 . The first group 
of authors claimed that as a result of the atomic vibra- 
tions in the /3-regime the boson-peak in a glass-system 
is observed, which depends on the fragility of the sys- 
tem. The more fragile the system ist, the weaker is the 
effects of the atomic vibrationen, i.e., the height of a 
boson-peak will be observed smaller and not clearer. As 
a result the relaxation-process is here observed clearer. 
In a connection with the fragility the glass-systems can 
be classified into two groups [32, |3a|: 1) fragile glass- 
systems, which have no directional bonds, e.g., Van der 
Walls and ionic systems. These glass formers exhibit 



strong non-Arrhenius-like increase of the viscosity upon 
colling; 2) strong glass-systems, which have strong direc- 
tional bonds. These glass-systems have the temperature 
dependence of the viscosity that is more Arrhenius like 
and weaker. 

In a connection with this classification it is until today 
still difficult to determine whether the metallic-glass be- 
longs to a strong glass-system or to a fragile glass-system, 
because usually the behavior of the viscosity in a metal- 
lic glass-system lies between both borderline case. So far 
from in hand analysis of our data we can state nothing 
about the "Boson-peak" . 

We observed that the self-part of the intermediate 
scatteing function for T < Tc shows a small bump at 
t w 2.5x10"^^ s. Th same phenomen war observed on 
MD-simulations for a OTP-system by Lewis and Wahn- 
strom [17] and for a binary Lenard- Jones mixture by 
Kob and Andersen [16]. This bump is interpreted by 
those authors as a result of the finite-size effect. The 
same bump was also observed on another Lenard- Jones 
mixture-system [SJ] , on a liquid-salt system [33] , and on 
a colloidal-suspenstion system [sg j. 

We have also analyzed that the height of plateaus as 
well as the time scale, at which the self-part of the inter- 
mediate scattering function decay finally to zero, depend 
strongly on wave vector q. Our MD-results of <&(q,i) for 
the longer time show that $(g, t) always decay to zero at 
both T <Tc and T > Tc. We interpret this result within 
the scope of the extended schematic MCT. The extended 
schematic MCT predicts that the long-time behavior of 
^{q,t) always decay to zero, when thermally activated 
hopping-processes is taken to be account in a system. 
According to Teichler [13, [H, [13] and Aspelmeier [U 
these hopping-processes take place on a MD-simulated 
NisoZrsQ. As one can there observe, the reason for dc- 
caies of ^{q,t) on the longer time scale is that on this 
time scale thermally activated hopping-processes run off. 
To prove this statement, we have to investigate the selft- 
diffusion konstant of atoms in our atoms. The MCT pre- 
dicts that self-diffusion konstant follows a power- law that 
is a function of Ty. 



Da oc {T ~ Tc 



(30) 



where 7 — l/(2a) + 1/(26) is valid and by knowing of 
the exponent A can be calculated, a is the specific 
atom. The estimation of 7 ove that expression is only 
possible, if a g-independence exponent A is given as re- 
quired by deriving that expression. On MD-simulation 
we can determine for different temperatures of interest 
the self-diffusion konstant throught the calculation of 
mean squared displacement (MSD). Then we find the 
curve Da vs T. The effective exponent 7 can be deter- 
mine by using the least-squared method. From the curve 
we can prove whether the thermally activated hopping- 
processes run off at T < Tc or tempertures near Tc in our 
system or not, as we investigate whether there is a de- 
viation from (|30p or not. If there is this deviation, then 
it means that the hopping processes still take place at 



T <TcOiT near T^ [4d, |43. 

From above discussion we can state that the behavior 
of ^{q,t) in our system agree with the prediction of the 
extended MCT, in which the long-time behavior of $(5, t) 
decays to zero as a resuh of the atomic diffusion. 

The MCT predicts that the behavior of $(g,t) in 
the last a-regime can be good approximated by the 
Kohlrausch-law. We have determined the Kohlrausch's 
parameters ^Oj Pi and tq. Our results show that the 
parameter in our system (see Table H]) depends weakly 
on wave vectors q and varies with temperature. The 
/9-value increases with a decreasing g-value about 0.05. 
This result is similar as one by Lewis and Wahnstrom in 
MD-simulated OTP-system ^^. 

On our results we have noticed that the /?- value seem to 
increase to one at higher temperatures with a decreasing 
g- value. This phenomane is also found in another MD- 
simulation by Bernu et.al. [43|. The MCT predicts that 
/3-value for <? ^ go (~ 20 nm~^), which corresponds to 
a nearest-neighbor distance, is found typically in range 
0.6 < /3 < 0.8 that depends on the system of interest; for 
(7 -C go the /3-value seem to increase. Also our /3-value 
agrees well with the MCT"s predictions. 

According to MCT the a-relaxation time diverge near 
Tc after the power-law [H, \I^ 



TABLE L The Kohlrausch parameters Aq, tq, /3 fitted to 
^{q,t) from MD-data for wave vector gg and qa- 



To{T) ^ {T - T.y 



(31) 



where 7 is the same exponent in Eq. (|30p . On our results 
we have found that the behavior of the power-law was 
broken in the nearest environments from Tc, because as 
a result of thermally activated atomic diffusion the a- 
relaxation time tq kept to a finite-value. 

Through a viewing of the a-relaxation time tq we have 
found that this relaxation time has a strong depending on 
g-value. The same result was observed in a MD-simulated 
binary Lenard-Jones mixture by Kob and Andersen |16| 
and in a experimentel quasy-elastic neutron scaterring 
on three polymers by Colmenero et.al. [4^. We have no- 
ticed that there is an anomaly on the a-relaxation time, 
namely, the relaxation time increases drastic near the 
glass-transition. One explains this anomaly [46] : the sys- 
tem is found above Tc on the way to a metastable equi- 
librium of undercooled liquid-phases. Then the system 
comes below Tc into a unstable state, where the relax- 
ation processes take place in a direction of the equilib- 
rium. On grounds of this anomaly we can state that the 
glass transition actually is a dynamic transition which is 
more in the change of the art on atomic movements and 
less in the change of the structure [4l|, |43| 

Figures S] and E] present the dynamic susceptibility. 
From both figures we found that there are three different 
frequency-regimes of the dynamics susceptibility as men- 
tioned in theory section. As predicted by extended MCT 
the a-peak shift with a decreasing temperature to the 
lower frequency, und this peak and the minimum of the 
dynamics susceptibility are no more dissapear at temper- 
ature near Tc- The extended MCT assumes this behavior, 
when thermally activited hopping-processes are included 
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~ 4.871 0.741 
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- 31.100 0.756 


0.911 
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formally in the memory kernel F{t). These processes has 
taken place in our system, i.e., these predictions agree 
with statements mentioned above. 

We have noticed that the height of a-peak for each 
temperature depends strongly on g-value, namely, the 
height of a-peak decreases with a increasing g-value. 
The microscopic (phonon)-regime as well as the height 
of this regime decreases with a decreasing both g-value 
and T. This result could be simply interpreted, if one 
takes into consideration the height of plateaus (the non- 
ergodicity parameter fc) on the intermediate scattering 
function, which depends strongly on the g-value. Because 
the height of the a-peak is proportional to the height 
of plateaus, and the microscopic-peak is proportional to 
1 — fc, then the above mentioned dependence of the both 
height-peaks results directly from the g-dependence non- 
ergodicity paramater. 
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FIG. 4: The dynamic susceptibility for wave vector qs 

Figures [6] and [7] present our evaluated kernels F{t) 
from MD-data of $(5, t) for both wave vectors of interest. 
From both figures we found that there is a temperature- 
dependence threshold-value $0(7") ( $o(r) £ [0; 1] ), 
where under ^o{T) the behavior of F{t) shows a course 
that corresponds to a polynom with positive coefficients, 
and otherwise this behavior deviats from a polynomial 
form. According to the prediction of the idealized MCT 
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FIG. 5; The same as Fig. |4l but for wave vector qg 



FIG. 6: The kernel F{t) for wave vector qs- 



the kernel F(i) as well as /($) should show a polynomial 
form. Also our results do not agree with this prediction 
for all of regimes, but only for the regime that lies un- 
der $0 {T) ■ This is understandable, because the idealized 
MCT does not fully describe the atomic vibrations of the 
system. 

We have noticed that ^o{T) depends strongly on the 
q- value and weakly on the change of temperature. ^q{T) 
shift with a decreasing both q- value and T to higher val- 
ues. 

There is here a important point that we can stretch, 
namely, by a calculating of the kernel F{t) we can make 
for $(T) < ^o{T) a new model of the memory kernel 
/($), which presents of course a polynomial form. By a 
extending for <i>(r) > ^o{T) we can then determine the 
term, which could be assigned to the atomic vibrations of 
the system and corresponds to the difference between the 
actually kernel F{t) and the fitted memory kernel /($). 
To this purpose we present our calculating for the kernel 
F{t) as a function $ in Fig. [5] and [D According to 
MCT one define /(<&) as follows : 



F<t) 



(32) 



/(a>) = ^A„$" 



where A„ are positive coefficients. We have fitted Eq. 
to the kernel F($). (see Fig. [TOl and [TT|) . and then we 
have for our system positive coefficients A„, which are 
given in Tables [TTl and UlTl 

To calculate Eq.®, as in I, we introduce the following 
equations 



G($) :=i^($)(I/$-l) 



.9($) = P($)(l/$-l) 



(33) 



(34) 



Eq.([34]) is analog to Eq.®, where /($) in g{^) is sub- 
stituted by P($). When we have calculated for $(r) < 
^oiT) F($) as well as /($), i.e., in compfiance with both 
equations we also have calculated for <I>(T) < ^o{T) the 
function P($). Therefore we can now state that P{^) 
is also a polynomial function with positive coefficients 
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FIG. 7; The same as Fig. [6] but for wave vector qa. 



A„, which will fit for $(r) < $o(r) the function g($) to 
G(<1>). 

In Figs. fT2l and [T3lwe present the G(<i?)-function. From 
both figures we can see that G(<I>) has a maximum at 
^o{T). This maximum shifts with a decreasing g- value 
to a higher value. We can now fit the function g($), 
which contains the polynomial function P($) with pos- 
itive coefficients (see Tables. HIl and Hill) , to G(<1>). Our 
results is shown in Figs. [14] and [151 From both figures 
we have for each temperature a maximum gm from (?($)- 
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FIG. 8; The kernel as a function of $ for wave vector qg. 




FIG. 9: The same as Fig. [8] but for wave vector ge- 
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FIG. 10: The kernel as function of $ for wave vector gg (dot- 
ted line: extrapolated low-<l> polynomial). 



function at $ = $„. This maximum is smaller als 1 for 
T > 1100 K ( g,n < 1 ), and for T < 1000 K nearly equal 

1 ( 9rn ~ 1). 

As mentioned in the theory section, gm describes a 
signature, so that, if gm < 1, the system is in a ergodic 
state (in a liquid state), and if gm > 1, it is in a non- 
ergodic state (in a structural arrest). Our results show 
that the Ni2oZr8o-system is for T > 1100 K in a liquid 
state and for T < 1000 K nearly in a structural arrest 



TABLE II: Coefficients of the polynomial expansion /($) 
'^^ A„<I>" for wave vector qs- 

An 



TK 


n 800 


900 1000 


1100 


1200 


1300 


1400 


1500 


1 0.255 


0.423 0.232 


0.357 


0.271 


0.253 


0.308 


0.198 


2 0.714 


0.794 0.658 


0.688 


0.661 


0.598 


0.820 


0.800 


3 0.040 


0.597 1.757 


1.177 


1.174 


1.630 


1.165 


1.592 


4 1.390 


1.318 2.144 


1.621 


2.089 


1.850 


1.876 


2.628 


5 2.218 


1.930 1.637 


1.655 


2.535 


2.116 


1.439 


0.057 


6 1.586 


1.628 0.700 


1.817 


2.414 


1.828 


0.593 


1.769 


7 0.321 


0.621 0.000 


1.445 


0.232 


0.566 


0.742 


0.040 


8 0.000 


0.000 0.009 


1.251 


0.012 


0.000 


2.089 


0.185 


9 2.438 


1.575 1.552 


0.623 


0.032 


0.000 


0.002 


0.143 


10 9.864 8.069 5.907 


2.253 


0.040 


0.020 


0.000 


0.451 



TABLE III: Coefficients of the polynomial expansion /($) 
X^n -^n^" for wave vector qa. 

An 



TK 


n 


800 


900 


1000 


1100 


1200 


1300 


1400 1500 


1 


0.235 


0.316 


0.362 


0.246 


0.173 


0.189 


0.240 0.162 


2 


0.021 


0.386 


0.373 


0.246 


0.394 


0.278 


0.394 0.135 


3 


0.918 


0.924 


0.990 


0.965 


0.214 


0.612 


0.682 0.515 


4 


1.291 


0.994 


1.240 


1.006 


0.529 


0.810 


0.640 1.397 


5 


0.016 


0.588 


0.873 


0.127 


1.448 


1.323 


1.563 2.250 


6 
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0.130 


0.289 


0.009 


2.594 


2.528 


2.250 4.494 
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2.873 


2.756 0.689 
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0.249 
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3.671 
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1.166 0.063 
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3.152 
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2.365 


0.020 


0.003 


0.003 0.001 


11 
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0.092 


0.030 
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1.323 0.003 


12 
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0.002 
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0.005 


0.001 0.014 
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FIG. 11: The same as Fig. [TD] but for wave vector qa. 



FIG. 12: G($) for wave vector qs according to Ea. (|33 



Ci.7 

O .5 

CO 



eCMD K 
SCMD K 
-1 OOO K 
1 lOO K 
1 £00 K 
1 300 K. 

1 500 K 




(1025 ± 25) K. 



FIG. 13; The same as Fig. [12] but for wave vector qa. 
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FIG. 14: Fitting 5($) (dotted lines) to G(<1>) (full lines) for 
wave vector qs according to Eqs.(|34[) and ((33}. 



with a correlation decay, which is as a result of a ther- 
mally activitated atomic diffusion. In Fig. [16] we present 
a curve gm vs T for each g-value. From that Fig. one 
can see that our system has a critical temperature Tc be- 
tween 1000 K and 1100 K. Here we decide Tr as about 



FIG. 16: gm as a function of T. 



B. Gaussian Approximation 

In following Figs. [H] [H [HJ and [ID] we have cal- 
culated the self-part intermediate scattering function 
^{q,t) by using Gaussian approximation, also according 
to Eqs. pTj) and p^ . As input data for both equations 
are the mean squared displacement (MSD) and the ve- 
locity autocorrelations functions obtained from our MD- 
Data. From figures we found that the gaussian approxi- 
mation is correct in the short-time regime, here also for 
t < 0.2 ps. For t > 0.2 ps the course of $(g,t) resulting 
from this approximation deviates from both MD-results 
and MCT- analysis- results. As mentioned in Subsection 
IIIIBI we have assumed that the DOS of phonons is equally 
approximated to the spectral distribution of velocity au- 
tocorrelations function (VACF). With this assumption 
we can state now that the phonon-regime in our system 
is the regime taken place in a range time smaller then 0.2 
ps. This time is a nearly temperature-independent time. 



Ct.9 




1 jC 



FIG. 15: Fitting g($) (dotted line) to G(<I') (full lines) for 
wave vector qg according to Egs.dM]) and ([55]) . 
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FIG. 17: The self-part of the intermediate scattering function 
for T = 900 K, according to Eqs. ([T7]l(line) and ([18} (dotted 
line), MD-result (dotted - dashed line), MCT-analysis (dashed 
line) . 
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FIG. 18: The same as Fig[T7]but for T = 1000 K. 
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FIG. 19: The same as Fig[TI]but for T = 1100 K. 



To gain more our statement that this regime is as a 
result of the phonon vibrations we borrow a important 
relation of the liquid theory, namely, the relation between 
the spectral distribution of VACF and the susceptibility 
of the self-part intermediate scattering function <i>(q,i). 
According to theory that relation is app roximated for a 
small g- value (lim (7 ^ 0) as follows |25l . l2a | 



x"{uj)^uj<^,{uj)^q^Z{uj)/uj 



(35) 



Fig[2T] presents our calculating results of Eg. ((35)) with 
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FIG. 20: The same as Fig[T7]but for T = 1200 K. 



an assumption that our gg-value is small enough. The 
calculating results show a fast nearly good agreement 
with that from the MCT-analysis (or a directly Fourier- 
tranformed $(g,t)). 
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FIG. 21: The susceptibility obtained from the spectral distri- 
bution of VACF through Eq.JSSI (lines) and from the MCT- 
analysis (dashed lines) 
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FIG. 22: The susceptibility obtained before superposition for 
wave vector gg. 
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FIG. 23: The susceptibility obtained after superposition. 

Actually we can do a inverse Fourier-transformation 
of the susceptibility obtained from Eq. ([55)1 to find back 
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FIG. 24: The master curve of the susceptibihty for qs ac- 
cording to Eq. (P^ . 
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FIG. 25: The new susceptibility after scaling for qs according 
toEa.(l27ll. 



FIG. 27: The susceptibility obtained after superposition. 



Ch.cr9 



r curvie 







FIG. 28: The master curve of the susceptibility for qe accord- 
ing to Ea. (f28)) . 



the self-part intermediate scattering function $(g,t), but 
from the result shown in FiglUwe do not do it, because 
Eq. (j35p is a approximated relation that is really correct if 
the wave vector q is very small (lim g ^ or, as one calls 
it, g is in a hydrodynamics limit), and then if we attempt 
to back Fourier transform of its susceptibility, surely we 
can not find a better <i>(g, t) as it found through Eq. pS)) . 



Further discussion about susceptibility we attempt to 
separate the susceptibility of this phonon-regime from the 
full susceptibility. As mentioned in Sec. |TTT]we assume 
that this separation is really to find a purely susceptibility 
of this regime. Our result for two wave vector is shown in 
Figs. [21] and [221 We scale now these susceptibilities, the 
result is shown in Figs [23] and [27l here we have chosen as 
a arbitrary normalization of temperature, Tq, 800 K, and 
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FIG. 26: The susceptibility obtained before superposition for 
wave vector q^. 



FIG. 29: The new susceptibility after scaling for ge according 
to Eg. pTl) . 
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we average these scaled susceptibilities through Eq.(3.14) 
to obtain the master curve of the susceptibility in this 
regime(see Figsi24land [28| . By using Eg. ([27]) we obtain 
the new susceptibility of this regime (see Figs. [55] and 



Adding the new susceptibility as shown in Figs [25] and 
1291 with the susceptibilities from other regimes, then we 
obtain the full susceptibility in the short-time regime as 
shown in Figs. [30] and [31] From both figures we found 
that the the peaks of the new full susceptibility for tem- 
peratures 1000 K and 1100 K are higher then those old 
ones, but that peaks for temperatures 800 K and 900 K 
are lower then those old ones. 
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FIG. 30: The full susceptibility before (dashed lines) and 
after scaling(full lines) in a short-time regime for wave vector 
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FIG. 32: The self-part of the intermediate scattering function 
in a short-time regime, MD-resuIts (Symbol), before (dotted 
line) and after scaling (dashed- line) for wave vector qs- 
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FIG. 31: The full susceptibility before (dashed lines) and 
after scaling(full lines) in a short-time regime for wave vector 



Taking back Fourier transformation of <i>c('j-') of the 
new full susceptibility in the short-time regime, then we 
obtain the new $(g,t) as shown in Figs. [5^ and [551 The 
deviation of the new $(g,t) from the old one for each 
temperature of interest is taken place at a time range 
that is smaller then 0.5 ps. 

It is interesting to calculate the new kernel F{t) and 
the new G(<l>)-function in this short-time regime. Our 
results for the kernel F{t) are shown in Figs. [Ml and [551 



FIG. 33: The self-part of the intermediate scattering function 
in a short-time regime, MD-results (Symbol), before (dotted 
line) and after scaling (dashed-line) for wave vector q^. 



The new kernel F(t) for each temperatures of interest 
show more regular then the old one. The deviation of 
the new kernel from the old one is taken place at a time 
range that is smaller then 0.5 ps, also this time range is 
the same as it found in the case of $(g,i). It is clear 
that $((?, t) in this regime depends strongly on the kernel 
Fit). 



V. CONCLUDING REMARKS 

Our system behaves as predicted by the schematic 
MCT-Modell, in the sense that the self-part intermediate 
scattering function $(5, i) shows at lower temperature 
three step of the structural relaxation. The behavior 
of the self-part intermediate scattering function $(g,t) 
agrees well for the long time with the prediction of the 
extended schematic MCT-Modcl in the sense that the 
self-part intermediate scattering function $(g,i) always 
decays to zero after long times as a result of thermally 
activated atomic diffusion. 

The behavior of the memory kernel F{t) agrees with 
the prediction of the extended MCT-Model in the sense 
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FIG. 34: The kernel F{t) in the short-time regime before 
(full lines) and after scaling (dotted lines) for wave vector qs- 
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FIG. 35: The kernel F{t) in the short-time regime before 
(full lines) and after scaling (dotted lines) for wave vector qe- 



that the behavior of F{t) below $o(i) corresponds to a 
polynomial function that has non-negative coefficients. 
This behavior depends strongly on the wave vector q. 



As in I, to determine approximately the critical tem- 
perature Tc we have used the maximum g„i of the char- 
acteristic function g($). We obtain Tc « (1025 ± 25) K 
for our system. 



Through the temperature superposition method in a 
short times regimes our results show that the vibration- 
phonon regime is really as a result of the harmonic 
phonon approximation, and then this regime can be good 
described by a gaussian approximation. In the case of 
a gaussian approximation, which follows from the liq- 
uid state theory, the courses of the self-part intermedi- 
ate scattering function ^{q,t) resulted from both VACF 
and MSD shows at time, t « 0.5 ps a deviation from 
the proper courses obtained from MD-simulations and 
MCT-analysis 
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